Energy of the nearest neighbor RVB state by systematic loop expansion 
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We present an approximation scheme for the calculation of the norm and energy of the nearest- 
neighbor-RVB state for the Heisenberg antiferromagnet on the 2D square lattice. The approximation 
leads to a systematic expansion of norm and energy, with the 'expansion parameter' being the 
maximum length of loops taken into account in the calculation of energy and norm. The expansion 
converges well, the best estimate for the ground state energy/site is — 0.434473J. 



The 'RVB spin liquid' is a frequently occuring phrase in 
connection with cuprate superconductors. It is supposed 
to describe a state with strong short range singlet correla- 
tions, but no long range magnetic order whatsoever. De- 
spite its frequently being referred to in the literature, up 
to now the 'RVB spin liquid' is a rather elusive concept. 
For example, there exists to date no simple and manage- 
able trial wavefunction for the 2D Heisenberg Antiferro- 
magnet which would describe such a state, nor is the pre- 
cise nature of its low lying elementary excitations known 
to any degree of certainty (as opposed e.g. to the antifer- 
romagnetic state, where linear spin-wave theory gives an 
excellent approximation to the ground state wave func- 
tion and allows a quantitative discussion of its excitation 
spectrum). Perhaps the best-defined 'RVB spin liquid' 
is the nearest neighbor RVB state - at least this 

wave function can be written down in compact form. We 
introduce the singlet generation operator on the bond i, j 
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where cj „—c\ CT Cj g.c\ ff , and the operator 



i+yj 



where i + x denotes the nearest neighbor of sites i in 
x-direction. The nearest neighbor RVB state on a 2D 
square lattice with N sites then is defined as 
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In other words, the state consists of a superposition of 
all possible distributions of N/2 singlets over the plane, 
<p, with equal phase. The total number of different sin- 
glet configurations corresponds to the number of dimcr 
coverings of the plane, N c . The calculation of N c is a 
well-known exactly solvable problem from statistical me- 
chanics [^J^] - however, knowledge of this number is not 
really necessary to estimate the ground state energy. As 
will be seen below, the energy of the state (|lj) is signif- 
icantly higher than that of the antiferromagnetically or- 
dered ground state - the nearest-neighbor RVB state itself 



therefore is no candidate for the ground state wave func- 
tion. On the other hand, and this is the main motivation 
for the present work, it has recently been demonstrated 
for Heisenberg ladders || that by starting from a 'sin- 
glet vacuum' a theory for triplet-like fluctuations 'on top 
of this vacuum can be constructed, which in fact com- 
pares very well to numerical results. While for ladders 
the topology of the system rather uniquely determines 
one specific singlet covering, the nearest-neighbor RVB 
state would be an appealing candidate for a generaliza- 
tion of this approach to an isotropic and translationally 
invariant state of a planar system. In trying to do so, 
however, one encounters an as yet prohibitive technical 
obstacle: the fact that singlet configurations correspond- 
ing to different dimer coverings are not orthogonal [0. 
More precisely, whenever we can draw a closed loop which 
passes through both sites of each dimer it touches, it is 
possible to shift the 'train of dimers' along the loop by 
one lattice spacing, so as to obtain a different covering. 
Then, the 'Neel components' in the two dimer coverings 
along the loop are identical, so that there is a nonvan- 
ishing overlap. Each 'Neel component' has the prefac- 
tor 2~'/ 2 , where I is the number of sites in the loop, so 
that the overlap originating from the loop is 2~('/ 2 ~ 1 ) 
B . Clearly the most important contributions thus come 
from states which differ only by short loops. First, the 
overlap due to short loops is larger, and second, a long 
loop 'blocks' more possibilities for drawing other loops. 
In the following, we want to make these considerations 
more quantitative. 
To that end we write the norm as 
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where ~ denotes the statistical average over all configura- 
tions \ip). As a first step we restrict the sum over <j) to 
those configurations which can be obtained from by 
'rotating' a certain number of ?=4-loops. We denote the 
number of these 4-site loops by v, and start with the case 
v=\. If there are N sites in the system, we can obviously 
draw N different 4-site loops in the plane. However, in 
most cases the dimers of \4>) will not 'fit' into the drawn 
loop. Let us assume that we have drawn a 4-site loop 
somewhere in the plane, and consider the lower left cor- 
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ner of the loop (see the square labelled a in Figure [l]) . 
The probability that the dimer which covers this lattice 
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FIG. 1. Interplay of loops and dimer coverings. 

site points into positive x or positive y direction, and thus 
'fits' into the loop, is 2/ z. Next, let us assume that the 
dimer on the lower left corner points in positive x direc- 
tion, and consider the upper left corner (see the square 
labelled b in Figure [l]. The dimer which covers this cor- 
ner cannot point into negative y direction, because the 
lower left corner is already occupied, ft thus has only 
z — 1 possibilities for its direction. If we now assume, 
that these z — 1 orientations have equal probability, the 
total probability that the dimers in the state l^) fit into 
the prescribed 4-site loop is z ^-\) ■ Whereas for peri- 
odic boundary conditions the first factor of 2/z follows 
rigorously from symmetry, the second factor of l/(z — 1) 
is an approximation. Numerical calculations on 4 x 4 
and 6t x 6 lattices with periodic boundary conditions 
indicate that the probability is higher, i.e. one should 
replace l/{z— 1)=0. 33333^0. 47. On the other hand the 
values of N c obtained for these systems differ strongly 
from the ones expected in the thermodynamical limit, so 
that the relevance for the infinite system may be ques- 
tionable. Since the difference is not too large anyway, we 
stick to our simple estimate. As will be seen below, a 
change of the probability can be incorporated easily. 
Next, if the dimers fit into the loop we can create a new 
state \4>) by rotating the dimers along the drawn loop, 
and have (0|^>)=l/2. The total contribution to the over- 
lap from all states which can be obtained by rotating one 
4-site loop in therefore is AAi with 
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Next, let us assume that we have drawn v different 4-site 
loops in the plane. The first question then is: in how 
many different ways can we do that? There are N ways 
to draw the first loop and since the second loop should 
not intersect the first one, a number of positions for the 
second loop are blocked. More precisely, the number of 
ways to draw the second loop is only N — bn, where 
&n = 9. If the first two loops are far from each other 
they will block 2bn possible positions, so that the num- 
ber of ways to draw the third loop is N~2bn. Continuing 
we estimate the number n„ of different ways to draw v 



loops as 
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The factor of 1/V! is due to the fact that by drawing one 
loop after the other, each configuration of loops is gen- 
erated in v\ different ways. We have assumed that each 
loop blocks precisely b\\ following loops, which is only an 
approximation; deviations will occur if a loop is placed 
'close' to a previously drawn loop. The probability for 
this to happen is proportional to the 'average density' of 
loops squared, and we will have to check later on that 
this is small. 

Assuming that (|J) is a reasonable approximation for the 
number of different loop configurations, we can now re- 
peat the considerations for the single loop for each of the 
v loops. In this way we find that the total overlap of \ip) 
with all states which differ by rotating v 4-site loops is 
simply n v X\, whence the contribution to the overlap of 
\i\>) with all states which differ by rotating an arbitrary 
number of 4-site loops is 
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The extra 1 on the r.h.s. is negligible and will be dropped 
henceforth. Note that we have extended the upper limit 
of the sum up to N/bn, where our above estimate cer- 
tainly is invalid. As mentioned above, our derivation will 
be valid only if the 'average density' of loops is small, 
and our next objective is to check this assumption. To 
be more precise, we want to find out which order v in the 
expansion (||) gives the dominant contribution. Approxi- 
mating the factorials by Stirling's formula and nominally 
differentiating with respect to the number of loops, v, we 
find that the maximum contribution to the sum comes 
from terms with 
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Moreover, in the neighborhood of vq we find n v \\ cx 
exp(—N(v — ^o) 2 )- On the average, each drawn loop 
thus has 20 'free squares' around it, and terms which 
significantly different loop density give a negligible con- 
tribution. This implies first that our approximation of 
neglecting correlations between the loops is probably rea- 
sonably justified for the terms which give the dominant 
contribution to the norm, and second that extending the 
sum in ([!]) up to N/bn is justified as well. 
Having obtained a first estimate for the norm of the RVB 
state we proceed to compute the expectation value of H 
to the same approximation. To that end we introduce 
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the projection operator V — \ijj)(ip\ and Q = 1 
find 

(</>\H\il>) = {4>\{V + Q)H\^) 

= (<j>\ip)E + (<p\QH\ip) 
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where E denotes the expectation value of the Hamilto- 
nian in the state \ip). We now write the expectation value 
of the Hamiltonian as 
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and again restrict the summation to those states which 
can be obtained from \tp) by rotating an arbitrary num- 
ber of 4-site loops. The first term on the r.h.s. in (Q) then 
simply gives the energy of one dimer covering, — 3-/VJ/8, 
multiplied by the approximate norm. The second term is 
more involved and represents the 'true off-diagonal' con- 
tribution. Let us assume that we act with the exchange 
term along a bond (i, j) connecting two sites covered by 
two different singlets in \ip). Then one has 
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where t\ v a creates the a-component of the 3- vector of 
triplets on the bond (i,i'). Next, a singlet and a triplet 
have opposite time reversal parity, and a necessary con- 
dition for a nonvanishing overlap along one loop is that 
both states have equal time reversal parity 'along the 
loop'. We can conclude that only those states \<j>) can 
have a nonvanishing matrix element with \ip) which dif- 
fer from \tp) along a loop covering both bonds, («,«') and 
Let us assume that we have a dimer configuration 
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FIG. 2. A contribution to h! = (<j>\Q,H\ijj): starting two sin- 
glets in ^-direction (left state) and acting with the transverse 
Heisenberg exchange along the indicated bonds one obtains 
two states which are orthogonal to the initial state, but corre- 
spond to singlets in y-direction. Each of these channels gives 
a total contribution of — J/8 to h' . 

which fits onto a 4-site loop, and consider the matrix el- 
ement of H between this configuration and the 'rotated' 
dimer covering. Looking at Figure ^| one can read off that 
the total contribution of this loop to the matrix element 
is 



(4\QH\4>) = -2-2-Z, 

where the additional factor of 2 comes from the 'spin re- 
versed' processes. Recalling that the same configuration 
would have contributed a factor of 1/2 to the overlap of 
the two states we find the total off-diagonal contribution 
is given by 
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This is just the expectation value (<f>\ QH\ip) / (<fi\ip) for a 
single 4-site loop times the 'most probable number' of 
loops, compare (^) . Dividing by the norm we obtain our 
first estimate for the expectation value of the energy: 
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] = -0.422619iVJ. (8) 



Quite obviously, the correction due to off-diagonal pro- 
cesses to the energy is rather small. 
We proceed to the second step. Whereas we restricted 
the sums in the norm and expectation value of H to con- 
figurations which differed only by 4-site loops from 
we now enlarge this subset to comprise all states differing 
by 4 or 6-site loops. Let us assume that we have drawn 
v 4-site loops in the plane, and ask how many ways are 
there to draw v\ additional 6-site loops. Putting the first 
6-site loop we have to make sure that it does not inter- 
sect any of the 4-site loops and assuming that each of the 
4-site loops blocks a total of 612 = 15 possible positions 
of 6-site loops, we find that the first 6-site loop can be 
drawn in Ni=N — b\iv ways. Now, we proceed precisely 
as for the 4-site loops: every 6-site loops that we draw 
blocks bi2 — 15.5 other possible 6-site loops (a 6-site loop 
in x-direction blocks 15 6-site loops in x-direction and 16 
6-site loops in y-direction) , so that we estimate the num- 
ber of ways in which the v\ 6-site loops can be drawn 
as 
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Thereby the extra factor of 2 Vl is due to the fact that 
there are two different types of 6-site loops (in x-direction 
and j/-direction). By analogous considerations as for the 
4-site loop, we now estimate that each drawn 6-site loop 
gives a contribution of 
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to the overlap. Proceeding as above we find that the to- 
tal contribution to the norm from all states containing v 
4-site loops and and arbitrary number of 6-site loops is 
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Comparing this to the previous expression, n^X", we note 
that the inclusion of the 6-site loops has two effects: a) an 
extra prefactor of (1 + 2A2&22) JV ^ b22 which is independent 
of v and which will simply enter as an additional multi- 
plicative prefactor to the total norm, and b) a renormal- 
ization of Ai: 
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Neither modifications poses any problem within our for- 
malism and summing over v we obtain an improved esti- 
mate for the norm: 

N 6 = N C (1 + 2A 2 &2 2 ) A7b22 (1 + \ib n ) N/bl K 

To compute the energy to the same approximation, we 
note that the contribution of a single 6-site loop to the 
off-diagonal matrix element (<f>\QH\ip) is — 2 J [|J. Wo 
thus find the total off-diagonal energy of all states with 
v 4-site loops and an arbitrary number of 6-site loops is 
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Again, this replacement has two-fold effect: there is a 
^-independent term of — (47VJA 2 )/(1 + 2X 2 b 22 ) times the 
total overlap, which represents the gain in energy due to 
inclusion of the 6-site loops and enters as an overall shift; 
second, there is a term cx v, which we can be reabsorbed 
into a renormalization of J in the off-diagonal energy due 
to 4-site loops: 
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This takes into account that drawing a 4-site loop the 
number of possible 6-site loops is reduced resulting in a 
reduction of the energy gain from a single 4-site loop. 
Summing again over v and dividing by the norm we ob- 
tain an improved second estimate for the energy: 
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= -0.4353387V J. 
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The 'most probable number' of 4-site loops drops by 
w 15% from 7V/21=0.0476iV to 0.04037V, whereas the 
most probable number of 6 loops is 0.0194iV. Enlarging 
the subset of states used for the computation of energy 
and norm does have only a moderate influence, i.e. our 
procedure is reasonably convergent. 
In addition the procedure has an induction-like charac- 
ter: enlarging the sums in (|^) and ^ the subset of states 
leads to additional prefactors in the norm, a renormal- 
ization of the preceding A's, additional terms in the ex- 
pectation values of H and a renormalization of the ex- 
change integrals in the preceding energy contributions. 
Repeating everything including also 8-site loops gives an 
estimate for the energy of E = — 0.4344737V J, the num- 
ber of 4-loops drops by only 8%. The convergence thus 
is quite fast, with the difference between subsequent es- 
timates dropping by an order of magnitude in each step 
(actually the dependence is not monotonous in that the 
estimate including 8-loops is higher than the the one from 
6-site loops). We can conclude that taking into account 
only 4-site loops may be sufficient to get reasonably accu- 
rate results with little effort - this may be an important 
simplification in further applications of the state ([!]) . 
Comparing with results of numerical calculations [QJ|] 
which give an energy of f=a — 0.6J/site the agreement is 
obviously not too good. This is probably due to thefact 
that our simpel estimates for the probability to find a 
given loop in an arbitrary dimer configuration are not 
too accurate. 
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